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Abstract. In this note we present algorithms for computing Euclidean min- 
ima of cubic number fields; in particular, we were able to find all norm- 
Euclidean cubic number fields with discriminants —999 < d < 10''. 



1. Introduction 

This article deals with the problem of determining whether a given cubic number 
field is Euclidean with respect to the absolute value of the norm. The correspond- 
ing problem for quadratic number fields was solved in 1952, when Barnes and 
Swinnerton-Dyer showed (after much work done by various authors) that the fol- 
lowing list of discriminants of norm-Euclidean quadratic number fields is complete: 

d = -11,-8,-7, -4, -3, 5, 8, 12, 13, 17, 21, 24, 28, 29, 33, 37, 41, 44, 57, 73, 76. 

In the cubic case, Davenport proved that the number of norm-Euclidean complex 
cubic number fields (i.e. cubic fields with unit rank 1) is finite, whereas Heilbronn 
conjectured that there are infinitely many totally real cubic fields which are norm- 
Euclidean. We hope that the methods presented in this paper will eventually lead 
to a complete list of norm-Euclidean complex cubic fields, and that extended com- 
putations for real cubic fields will show whether Heilbronn's conjecture is reasonable 
or not. 

2. Notation 

In order to present our method we need a few definitions. Let be a number 
field, and let Ok denote its ring of integers. The Euclidean minimum of ^ G is 
defined to be 

MiK,0 = inf {\Nk/q{^ Ok}- 
The field K is Euclidean with respect to the absolute value of the norm (norm- 
Euclidean for short) if M(K,^) < 1 for all ^ G K. Let us introduce the Euclidean 
minimum M{K) of K by putting M{K) = sup {M{K,^) : ^ e K}. Obviously, K 
is norm-Euclidean if M{K) < 1, and not norm-Euclidean if M{K) > 1 or if there 
isa^e K such that M{K) = M{K, ^) = 1. 

Next we introduce the inhomogeneous minimum. To this end let K be a number 
field generated by a root a of an irreducible monic polynomial / G Z[a;]. Let 
ai, . . . , ar denote the real roots, and Or+i, ar+i, . . . , as, cEJ the s pairs of complex 
conjugate roots of / in C; then the maps a — >■ aj can be extended to yield r 
embeddings : K — > M and s pairs of complex conjugate embeddings 

<^r+l, 4>r+l, ■ ■ ■ ,(l>s,4's ■ K > C. 

Choose a Q-basis {f3i, . . . , /3„} of K; the map 



^ R" : 2^ ajPj ^ (ai, . . . , a„) 
1 
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embeds K into R", and we will identify K and 'k{K) for the rest of this article. 
Clearly K is dense in R" , so we will write K = M" if we want to make it clear that 
we regard M" as the closure of K. Now 

r s r» 

N:R"^R:{xu...,Xn)^l[\j2^iM^i)\- H 

j—l i j—r+1 i 

is a continuous map which coincides with the absolute value of the norm Nj^/q 
when restricted to K. Similarly, the maps 

I • |j : M" — ^ M : a; = {xi, . . . ,a;„) ^ \x\j = | , (0 < i < r + s), 

i 

are continuous and their restrictions to K agree with the r + s archimedean val- 
uations of K. By an abuse of language, we will refer to the maps N and | • \j as 
the 'norm' and the 'valuations' of respectively, although is not a norm on K 
since N{x) = does not imply x = 0. Similarly, the | • |j are not valuations of K 
for the same reason. All we can say is 

Proposition 1. Let K he a number field, and assume that limj_>.oo = for all 
^ < j < r + s and a sequence of elements G K. Then Wmi^ao = 0. 

Proof. K is an n-dimcnsional Q-vector space, with a nondcgencratc bilinear form 
given by {$,,rf) = Trx/Q(^77), where Tr^/Q denotes the trace of if/Q. Choose a 
Q-basis {ai, . . . , a„} of K, and let . . . , denote the dual basis with respect 
to (•, •), i.e. the basis with the property Trx/Q(ai/3j) = 5ij (Kroncckcr's delta). 

Now assume that limj_,.oo = for all 1 < j < r + s, and let 5 > be given; 
then there exists an TV e N such that < 5 for alH > A'' and 1 < j < r + s. 

Write as = x^^'ai + . . . + Xn^'oin- Then \x^^^ \ = \ Trj<-/Q(/3fc^i)|. Since the trace 
is the sum of all conjugates of (ik£.ii applying the triangle inequality yields 

\xf\ < \Pk^^\l + . . . + \M^\r + 2\/3k^,\r+l + . . . + 2\Pk^i\r+s 

< S{\(3k\i + ... + 2\(3k\r+s)<SC, 

where C = \pk\i + ■ ■ ■ + 2|/3fc|r+s does not depend on i or the choice of S. Since we 
can make S as small as we please, we find that limj_>.oo x^^^ = for all 1 < A; < n, 
and this is equivalent to limi^oo 6=0. □ 

Obviously K is norm-Euclidean if and only if for every ^ £ K there is an a e Ok 
such that — a) < 1. Actually, all known examples of norm-Euclidean number 

fields satisfy the stronger condition that for every ^ G K there exists an a € Ok 
such that N{(, - a) < 1. We put 

M(K,C)=mi {N{C-f]):f]eOK}, 

and define the inhomogeneous minimum of K as M{K) = sup {M{K,^) : ^ G K}. 
Clearly M{K) < M(K); it is conjectured that M{K) = M(K) for all number fields, 
but so far this equality has been proved only for fields with unit rank < 1. 

We say that M(i?) is isolated \i M2{K) = sup {M(F^^) : ^ e :^\{C : M(F,C) = 
M{K)}} < M(K). In this case, we call M(K) = Mi(K) the first and M2(K) the 
second minimum of K. 

Remark 1. K is a ring. In fact, the product of two elements ^ = ^Oia' and 
r] = ^ ha^ of K has the form ^ Cia\ where Cj is a polynomial in Q[ai, . . . , 
Thus, the product ^rj can be given a meaning for real values of the coefficients a,, 6,. 
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Now put n = {K : Q), and choose an integral basis {/3i = 1, (32, ■ ■ ■ , Pn} and a 
real number A; > (for example k = 0.99). We start by dividing 

n 

F+^{^^Y1 I «i ^ 1/2], 02, ... , an e (-1/2, 1/2]} 

into smaller subcubes. Such a subcube S' is called k-covered (or simply covered if 
the reference to fc is clear) if we can find a 7 € Ox such that A^(^ — 7) < fc for all 
^ e S*; it is called uncovered if we cannot find such a 7 (even if there exists one). 
Finally, a point £^ € K is called k-exceptional if — 7) > fc for all 7 e Ok. 

Remark 2. Observe that is 'half a fundamental domain' in the sense that 
every ^ € K /Ok has a representative in F = U where = —F+ (and 
only one unless the representative lies on the boundary of F). It is clearly sufficient 
to consider F+ since N{—^) = N{^). For cyclic cubic fields K we could reduce F+ 
further by exploiting the fact that Ni^") = N{^) for aU a € Gal{K/Q). 

Remark 3. Occasionally it simplifies computations to use fundamental domains 
other than F; they will be denoted by F, and we will always assume that F has 
compact closure. As an example, take 

n 

F=[( = Y^ a,f3, I ai e [0, 1), 02, • ■ • , a„ e (-1/2, 1/2]}. 
1=1 

Remark 4. For real numbers k',k > such that fc' > fc it is clear that any 
fc-covered set is also fc'-covered, and that any fc'-exceptional point is fc-exceptional. 

3. The Algorithms 

In this section we will describe the five programs (Eu3_l,. . . , Eu3_5) which have 
allowed us to compute Euclidean minima M{K) for many cubic number fields K . 
Since we used floating point arithmetic to compute the M{K) (which are rational 
numbers, at least in each case we succeded in its computation), a few explanations 
are in order. Suppose that we want to show M{K) = c for a number field K. 
Then we choose fc < 0.99c (as a protection against rounding errors), and, using 
the programs Eu3_l - Eu3_3, we compute cubes Sj which contain all fc-exceptional 
points (we want to be sure that they contain every c-exceptional point). Then we 
exploit the action of the unit group Ek on these cubes to compute the possible 
exceptional points, and since this is done with integer arithmetic, we are able to 
get exact results. 

Experiments with e.g. fields whose minima are known from hand computations 
(choosing values of fc close to the minimum and using very small cube lengths t) 
have led us to trust our results. Moreover, all our results agree with those obtained 
before (e.g. by Smith [IT and Taylor [12 ) or during the writing of this article (by 
D. Clark [5] and R. Queme). 

Our programs require as input a file called disc (i.e. 985 for the field with dis- 
criminant d = 985; for complex fields we used the absolute value of the discriminant 
preceded by for example _199 for the field with discriminant d = —199). This 
file contains the following data: 

• disc if; 

• the coefficients of the irreducible monic polynomial /; 

• — for real fields the roots of the polynomial a, a' , a"; 
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— for complex fields the real root a, the real and the imaginary part of 
the complex root; 

• the coefficients with respect to the base {-,—,—} (where g is defined 
below) of a system of independent units; 

• the index g — {Ok ■ '^[o]) and in the case oi g ^ 1 also gx, gy, gz, where 
{1, a, 61 = ^ + + ^a^} form a Z-basis of Ok] 

• the value k > 0; 

• the edge length £ of the cubes; 

• the coordinates of the uncovered cubes. 

For the sake of simplicity we treat only cubes having the same size; thus a cube 
is uniquely determined by its leftmost corner and the edge length £. We therefore 
start with the four cubes making up and the initial file 985, for example, looks 
as in Table 1 below: 





Table 1 




Table 2 




Table 3 




985 


1 


-6 


-1 


985 1 


-6 


-1 


985 1 


-6 


-1 


-2. 


93080160017276 


-2.93080160017276 








-0. 


16296185677753 


-0.16296185677753 




0.02 




2. 


09376345695029 


2.09376345695029 


0.38 


-0.22 


0.38 







1 








1 





0.38 


-0.22 


0.4 




2 


-1 





2 


-1 





0.38 


-0.2 


0.38 






1 






1 




0.38 


-0.18 


0.38 






0.9 






0.9 




0.38 


-0.2 


0.4 






0.5 






0.1 




0.4 


-0.22 


0.38 







-0.5 


-0.5 


0.3 


-0.5 


-0.5 


0.4 


-0.22 


0.4 







-0.5 





0.3 


-0.5 


-0.3 


0.4 


-0.2 


0.38 










-0.5 








0.4 


-0.18 


0.38 













0.4 


0.4 


0.4 


0.4 


-0.2 


0.4 



We now run the programs Eu3_l, Eu3_2 and Eu3_3, which will be described in the 
sequel, on the file disc. The programs Eu3_4 and Eu3_5 will be explained before 
Prop. Eland Cor. U\ respectively. 

Eu3_l. This program first asks for a discriminant, then reads the corresponding file 
disc. The first eight lines of disc are copied to the (temporary) file disc .new. The 
next input is an integer / which is the factor by which we divide the edge length of 
the cubes. We have used / e {1, 2, 4, 5}, depending on the size of disc; of course 
the choice / = 1 is only useful after k has been replaced by some k' > k. Thus 
Eu3_l reads £ from disc and writes £/ f to disc. new. Moreover, if no file disc.p 
exists, Eu3_3 creates one and writes the translation vector (0, 0, 0) into it. 

Now Eu3_l splits each cube read from disc into smaller ones, computes an 
upper bound B — B{S) of the minimum for each of these subcubes, and writes 
those S with B{S) > k to the file disc. new. Having reached the end of the file 
disc, it copies disc to disc.bak (a security backup) and disc. new to disc. 

How do we bound the minimum on a cube S = [a, a+£] x [6, b+£] x [c, c+i?]? Since 
the norm is the product of the three if -valuations, we only need to find bounds for 
I^Ij, where ^ = x + ya + z9 S. But since is a linear function of x,y,z, it 
takes its maximum at the corners of the subcubes. Instead of computing | • \ j at all 
eight corners and taking the maximum of these values as our upper bound (as the 
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programs with which we computed the tables at the end did), we can use a trick 
due to Roland Queme [10] which gives this bound at one stroke (but which doesn't 
seem to work except for valuations corresponding to real enibeddings): in fact, the 
value of I • Ij at the eight corners is one of 

|^o±|±|a±|0|j, where = (a + |) + + §)« + (c + 

the corners corresponding to the different choices of the signs. Using the triangle 
inequality we easily get 

|eo±|±|«±|e|, < ieol, + |(l + l«l, + l%). 

On the other hand, choosing the signs in such a way that £,0, ±1, ±a and ±9 all 
have the same sign (in the embedding corresponding to | • j^ ; here is where we need 
that the embedding is real) we see that this bound is best possible. 

Using this method of bounding N{S — 'y), Eu3_l will start looking for translation 
vectors in disc.p; if there is no element in this file such that N{S — 7) < k we 
search for translation vectors in the set 

(1) I = {x + ya + z0\{x,y,z) e I.\\x\ < M,,\y\ < My,\z\ < MJ, 

where Mx, My, Mz were usually chosen (depending on £) as follows: 





£ > 0.02 


0.01 >£> 0.001 


£ < 0.0005 




8 


19 


30 


My 


5 


12 


17 


Mz 


2 


3 


5 



The nonsymmetric limits were suggested by experience. If we find a translation 
vector 7 e Ok such that N{S — 7) < fc we write 7 to the file disc .p. There are 
several reasons for proceeding like this: 

(1) If N{S — -f) < k, then 7 has a good chance of satisfying N{S' — 7) < fc for 
cubes S' in the vicinity of S. By searching disc.p first we actually save 
much CPU time. 

(2) If we find that we have to replace k by some k' < k, we have to redo 
the computations from start; we are, however, able to use the translation 
vectors found in the previous runs. In fact, our programs allow the option 
of searching for new vectors or just using those in disc.p. 

After the first run of Eu3_l with / = 5 (which took 0.63 seconds of CPU time on 
an RS 6000), the file 985 looks as in Table 2. It contains exactly 106 noncovered 
cubes. The file 985. p contains the following translation vectors: 

(0 0), (0 10), (0 2 1), (10 0), (-1 -4 -1), (0 3 -1), (-1 -3 2), 
(0 1 -1), (0 2 0), (3 -1), (0 2 -1), (-1 0), (-1 -4 2), (-2 0) 

Eu3_2. This program acts like Eu3_l with the difference that the original cube is 
written to disc. new as soon as one of its subcubes cannot be covered. In other 
words: Eu3_2 eliminates those cubes whose subcubes of length £/ f can all be cov- 
ered. This is convenient if we already have to deal with a lot of cubes and a further 
division done as by Eu3_l is likely to lead to an enormous number of smaller cubes. 
We usually run Eu3_2 before using Eu3_3 in order to save CPU time. 
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Eu3_3. Scanning through ah the integers of / (see (HJ) takes much time. We can 
avoid searching for exphcit translations if we proceed as fohows: we muhiply a cube 
T with a non-torsion unit e, translate the result back into the fundamental domain 
by subtracting /3 e Ok and look whether eT — /3 intersects one of the cubes not 
yet covered. The program does not really compare the oblique prism eT — (3 with 
the uncovered cubes but rather uses the smallest box S which contains T and has 
faces parallel to the coordinate planes (this will be improved in the next version 
of our programs; we also remark that - in order to avoid rounding errors - we do 
not compare eT — /? with T and —T but with a slightly larger cube obtained by 
adding (resp. subtracting) to (resp. from) the coordinates of the corners of 
T). Evidently we have to compare the box also with the 'opposite' cubes, i.e. the 
cubes multiplied by —1, since we only kept the 'bad' cubes of half the fundamental 
domain F^. If there is no intersection, the cube itself can be eliminated. 

The reason why this program is so successful is the following: suppose that S 
is a subcube such that eS — a (where a € Ok is the element needed to translate 
eS back into F) is covered by 7 S Ok, i-e. N{eS — a — ^) < k. This means of 
course that N{S — (3) <k for (3 = e^^{a + 7): but /? will usually have much larger 
coefficients than those scanned in Moreover, in general eS — a will not be 
covered by a single element 7 S Ok, which means that we would have to divide S 
into subcubes before we could cover it directly. 

A run of Eu3_l on the file in Table 2, again with / = 5, leaves only 27 subcubes 
uncovered. Running Eu3_3 twice on the file obtained we are left with 10 uncovered 
cubes (see Table 3). Running Eu3_2 on the file in Table 3 deletes (0.38 -0.18 
0.38) and (0.4 -0.18 0.38). Now we have covered F except for the set T = 
[0.38,0.42] X [-0.22,-0.18] x [0.38,0.42]. Muhiplying the corners F of T by the 
unit a we find 



p 


aP 


0.38- 


0.22a- 


h 0.38a^ 


0.38- 


0.34a - 


f 0.40a^ - 


- 3a 


~a' 


0.38- 


0.22a- 


h 0.42a2 


0.42- 


0.10a- 


f0.36a2 - 


- 3a 


-a2 


0.38- 


0.18a- 


h 0.38a2 


0.38- 


0.34a - 


f 0.44a2 - 


- 3a 


-a2 


0.38- 


0.18a- 


h 0.42a2 


0.42- 


0.10a- 


f 0.40a2 - 


- 3a 


-a2 


0.42- 


0.22a- 


h 0.38a2 


0.38- 


0.30a - 


f 0.40a2 - 


- 3a 


-a2 


0.42- 


0.22a- 


h 0.42a2 


0.42 - 


0.06a - 


f 0.36a2 - 


- 3a 


-a2 


0.42- 


0.18a- 


h 0.38a2 


0.38- 


0.30a - 


f 0.44a2 - 


- 3a 


-a2 


0.42- 


0.18a- 


h 0.42a2 


0.42 - 


0.06a - 


f 0.40a2 - 


- 3a 


-a2 



This shows that aT is contained in the set 

T' = [0.38, 0.42] X [-0.34, -0.06] x [0.36, 0.44] 3a - a^ 

(in general, this is a very crude estimate; the next version of our programs will take 
the actual shape of eT into account). Observe that T' — 3a -I- a^ does not intersect 
—T. If we were to keep dividing the uncovered cubes, this picture would not change: 
the length of the cubes would become smaller and smaller, and so would the size of 
our uncovered set T, but aT — 3a -I- a^ would always have points in common with 
T . This is the situation which is described in 

Proposition 2. Let K he a number field and e a non-torsion unit of Ek- Suppose 
that T C F (vf. Rem. has the following property: 
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there exists a unique (3 g Ok such that, for all ^ G T, the element 

eS^ — 13 lies in a k-covered region of F or again in T. 
Then every k-exceptional point G T satisfies I^Co — ^zjIj ~ Q for every K -valuation 
I • Ij such that \e\j > 1. //, moreover, \e\j ^ 1 for all K-valuations, then the sequence 
Co, Ci: ^2, ■ • ■ of k-exceptional points defined by the recursion ^i+i = e^i — /3 satisfies 
lim^—j-oo C^ — £^T ■ 

Proof. Suppose that ^ T is fc-exceptional. Since M{K,^) — M{K,e£, — (3), our 
assumption imphes that all the defined by ^i+i = e^i — /3 are fc-exceptional points 
in T. Now C — is the fixed point of the map ^ i — > — /3. Induction shows 
that £,i — C — £'(Co ~ C) for * > 0. In particular we see that 

(2) 16 -CI, = \e\)\^o-C\r 

Now we claim that there exists a constant C > such that j^i — C|j < C for all 
i > (C = C{j) may depend on j, but we can always choose C as the maximum 
of the (finitely many) C{j)). In fact, since ^i, ( G F, we see that their difference is 
an element of 2F C K. Since 2F has compact closure and | • \j is continuous, | • \j 
has a maximum C on the closure of 2F and hence is bounded on 2F. 

Assume that \e\j > 1; then the fact that the left hand side of Equ. Q is bounded 
implies that |Co — Cli =0- This in turn gives immediately I'Ci — Clj = 0- If, moreover, 
\e\j 7^ 1 for all j < r + s, then either \e\j > 1 and — (\j = 0, or \e\j < 1 and 
= lim^^oo — C\j- This implies lim^ = C by Prop. [1] □ 

In our case there are three embeddings; we have \a\j > 1 for j = 1,3, and 
\a\j < 1 for j — 2. From Prop. [5] we can deduce that every 0.9-exceptional point 
C e T satisfies \^ - C|i = |C - CIs = 0, where C = ^^5f^ = i(2 - a + 2a'^). In order 
to show that ^ is in fact the only 0.9-exceptional point ^ g T we have to prove 
1^ — ^12 = 0; this is done by using the inverse of the unit a: 

Theorem 3. Let K be a number field, T a compact subset of F, and let e G Ek 
be a non-torsion unit. Suppose that 

(1) there is a (3 d Ok such that, for all ^ G T, — /3 lies in a k-covered region 
of F or again in T ; 

(2) for all ^ T there is a ^ £ Ok such that e^-'^^ — 7 lies in a k-covered region 
of F or again in T; 

(3) \e\, forl<j <r + s. 

Then is the only possible k-exceptional point of T . 

Proof. Let ^ G T be a /c-exceptional point. If \e\j > 1 then Prop. [2] shows that 
IC — TTtIj — 0. The other iiT-valuations | • \j satisfy \e\j < 1 because of 3., and 
we see \s^^\j > 1. Since ^ is /c-exceptional, so is d = £^^S, — 7 G T. Now 
- ^ = e{e-^^ -7)-^ = ^-(/3 + £7); but ^ € T and ^ - (/3 + £7) G T imply 
that f3 + e-f ~ 0, i.e. 7 ~ —e^^/3. Therefore the 7 in 2. is uniquely determined, 
and we can apply Prop. [5] with e^^ and 7 instead of e and (3. This shows that 
any fc-exceptional point ^ G T satisfies |C — ^_7_-^ \j = for all if- valuations with 
\e\j < 1. But = = Thus \^- =0 for alll < j < r -I- s, and 

by Prop. [l]this implies that ^ = j^j. □ 

Remark 5. This theorem is attributed to Cassels in [1]. 
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Remark 6. For every number field K there exists a complete system of indepen- 
dent units Ei such that \ei\j ^ 1. This follows directly from Minkowski's proof of 
Dirichlet's unit theorem. 

Remark 7. If condition 1. of Thm. [3] is satisfied but 2. is not (for example if 
there is a second uncovered subset T' such that e~^T and T' have common points 
mod Ok), then T might contain irrational exceptional points (converging to ^zy) 
as in the last sentence of Prop. [21 

Remark 8. If eT intersects — T, try to apply Thm. |3]with e replaced by —e. 

Remark 9. Suppose that (^, ^, \) is one of the A;-exceptional points; in this case, 
there will be uncovered sets in all eight corners of the fundamental domain F . In 
order to apply Prop. [2]one has to choose F e.g. as in Remark[3l because this allows 
us to collect these uncovered cubes into one set T lying in the center of F. 

In our example, computations similar to those above with a^^ = —6 + a + 
instead of a show that a~^T is contained in 

T" = [0.26, 0.54] X [-0.24, -0.16] x [0.38, 0.42] + 3 - a. 

Now Thm. 13] shows that the only possible fc-exceptional point (fc — 0.9) in T is 

^ 2 1^ 2 2^ 

^ OL~l 5 5 5 

Eu3_4. This program does the necessary computations: it checks whether a cube 
S multiplied by a non-torsion unit e and translated back into the fundamental 
domain F intersects either S or —S] in both cases, the possible exceptional point 
is computed and written to the file disc .n. We can also search for orbits of length 
> 2 by replacing e by e" for some n > 2, provided I is small enough. Verifying 
the conditions of Thm. [3] are currently still done by hand; see, however. Rem. 
1121 The actual computation of the possible exceptional point is done using integer 
arithmetic; whenever the precision was insufficient (e.g. for disc if = —680 or —728) 
we used PARI. 

The next question is how to compute M{K^S^). This is done as follows: first we 
notice that M(K,^) = M(K,^ + a) for a G Ok, i-e. M(K,^) only depends on 
the coset ^ of ^ in 'K/Ok- Next we observe that M(K, ^) = M(K,e^) for 

units e e Ek. For ^ e K/Ok put Orbe(C) = {^"^ : n G Z} (this is the orbit of ^ 
under the action of e) and Orb(^) — {e£, : e € Ek}- Then the Euclidean minimum 
is constant on every orbit. 

Proposition 4. For number fields K with unit rank > 1, the following properties 
o/C e K/Ok are equivalent: 

i) Orb(^) is finite; 

ii) Orb£(^) is finite for some non-torsion unit e € Ek; 

iii) e e K/Ok- 

Proof. The implication i) ii) is obvious. Assume that e is a non-torsion unit 
such that Orb£(^) is finite. Then there exists an n e N such that e"^ = ^, and this 
implies that ^ — + Ok for some a e Ok, i-c. £, G K/Ok- 

Finally assume that ^ = f+OK for some a, /3 G Ok- If G Ek, then ^ = 0+Ok 
and the claim is trivial. Otherwise observe that multiplication by a unit e maps 
^ to some element of the form %- + Ok, where a' = ae mod /3. This shows that 
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#Orb(C) < #{0k/I30k) = N{13). In fact, if Ok has class number 1 and if 
(a,/3) = 1, then we clearly have # Orb < ^{Ok/I30kV. □ 

In this paper we will not deal with computing minima for ^ e K /Ok with infinite 
orbit (the last sentence of Prop. [2] gives a hint as to how such infinite orbits might 
arise), so we assume from now on that ^ e K/Ok- The basic idea how to compute 
M{K,£/) is due to Barnes and Swinnerton-Dyer (see [T], Thm. B, for the case of 
real quadratic number fields). 

Proposition 5. Let K — Q(q!) he a number field with unit group Ek- If, given 
£^ G K and a real number k > 0, there exists 7 S Ok such that N($^ ^ 7) < k, then 
there exists — X]j=o o.jO^'' £ ^ with the following properties: 

(1) C = £,j mod Ok for some £_j £ Orb(^); 

(2) \ai\ < (0 < i < n) for some constants /i^ > depending only on K; 

(3) N{0 < k. 

Since the number of elements of K satisfying 1. and 2. is finite, we can prove 
M{K,£^) > fc by simply computing the norms of all these elements. We will prove 
Prop. [5] only for cubic fields. 

Let K = Q{a) be a cubic number field; replace ^ by ^ — 7 and choose a unit 
e G Ek such that the conjugates of /3 = = a + ba + ca^ are small. Let us consider 
the following system of equations: 

/3 = a + ba + co? 
13' = a + ba' + ca'^ 
P" = a + ba" + ca"^ 

This system is linear in a, 6, c, and the square of its determinant is 



1 


a 


a" 


1 


a' 




1 


a" 





A = det \ 1 a' a' = disc(l,Q;,Q; ), 

V 1 a" a"^ J 

which is clearly ^ 0. In fact, we have A — g^ disciiT for some integer g called the 
index of a. Therefore we get, by Cramer's rule, 

y/Aa = Pa'a"{a" - a') + P'a"a{a- a") + I3"aa'{a' ~ a), 

VAc = I3{a"-a')+l3'{a~a") + l3"{a'~a). 

In order to compute bounds for a, 6, and c we have to find good bounds for the 
conjugates of /3 = ^e. 

We begin with the complex cubic case; as for real quadratic number fields, there 
is one fundamental unit 77. Replacing 77 by 77^^ if necessary we may assume that 
1 77 1 > 1. Now for every ci > there is a unit e — rf"- such that 

ci < M < ci • |?7|. 

Since |^'e'| = |f"£"| in the complex case, we get 
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where we have put k = N{£_) = iV(^e). We will choose ci in such a way that the 
resulting bounds on X ^ \/3{a" - a')\, Y = - a")\ and Z = |/3"(a' - a)| are 

equal. A little computation shows cf — \ a'~^a''\'^ ' ^^^^ ^^^^ yields 

x,Y,z < vmvwi 

Applying Lemma |6] below to x = X,y — Y and z — Z we find xyz — k\J\K\ and 





|c| < 

The obvious bound 

|6| < (X|a"+a'|+r|a"+a|+Z|a+a'|)/%/A < 

can be sharpened by applying Lemma [5] to x — X\oi 
z — Z\a + a'\, and the same goes for \a\. The actual bounds coming from Lemma [S] 
are computed by machine in each case because they depend on the size of |q!+q;'|, . . . 
etc. This concludes the proof in the complex cubic case. 

Remark 10. The bounds in the above proof are much better than those obtained 
by Cioffari [1] for the case of pure cubic fields. 

Lemma 6. Suppose that a;, y, z e M satisfy the inequalities < x < ci, < y < C2, 
< z < C3, and < xyz = k. Then 



X + y + z < max < Ci 

Proof. We want to find bounds for f{x,y) = x + y + —. Now / is positive in the 
domain under consideration 

D = \x,y > : X < ci,y < C2,xy > — |, 

its gradient vanishes only at x = y = and its Hesse matrix there is positive 
definite; this implies that / takes its maximum on the boundary. 

Assume for example that x ^ ci; then we have to find an upper bound for 
/i (y) = ci+y + . Again, /i assumes its maximum on the boundary. For y — C2 
we get the bound ci + C2 + ^r^; on the other hand from z ~ < C3 we get 
2/ > and we find /i(^) = ci + cg + 

The cases y — C2 and xy — ^ are treated similarly. □ 

Let be a real cubic number field, and let rji and 772 denote two independent 
units. We will denote the conjugates oi ^ £ K hy ^, ^' , ^" . For units e E Ek and a 
fixed embedding K — > R we define 

" kl, if |e| > 1, 
\e\~\ if|e|<l, 



lie) := 



and we put 71 = 7(771), . . . , 7^ = 7(?72)- 

Suppose that ^ £ K has norm k\ then, for any real numbers ci,C2 > 0, we can 
find a unit e G (771,772) such that (cf. ^9 ) 

Ci < |^£| < C17172, C2 < IC'e'l < £27^72- 
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This gives us the fohowing bounds on 

Now we proceed as in the complex case, put (3 = ^e, and choose ci, C2 in such a way 
that the resulting bounds on X = |/?(a"-Q;')|, y = and Z = |/3"(q;'-q;)| 

are equal. In fact, setting 

3 _ u - «'ll« - 3 _ 7i72|a - - a"| 

•^1 ~ ^ ,,2_,2|„,/ „,HI2 ' ^2 — K- 



yields the bounds X, F, Z < \/k ■ 7i727i72 VA. In particular, we find 



)={X + Y + Z) < 3\fk^^^i^^i^^i^^. 



"=1 ^ ./A 

Making use of Lemma|n]we can improve this by a factor of almost 3/2. The bounds 
for \b\ and \a\ are derived similarly; this concludes the proof of Prop. [S]in the cubic 
case. For general number fields, the proof makes use of the dual basis (cf. the proof 
of Prop. [J). 

Let us get back to our example of real cubic field with discriminant d — 985. 
Let a denote a root of /. Then {l,a,a^} is an integral basis of Ok, and two 
fundamental units are given by 771 = a and 772 = 2 — a. Put ^0 = -1(2 — a + 2q;^); 
then Orb(^o) = {Co}- Using k = 1.05 wc get the bounds ^1 — 6.2, ^2 = 3.2, 
/^3 = 1-9. We find M{K,^o) = N{^o - 2) = 1. 

Eu3_5 is the program which does these computations. In fact, for each ^ in our list 
disc.n, it calculates 

min {N{^j + a + ba + c6)} 

for all (a, b, c) S Z"^ such that the coefficients of the sum +a + ba + c9 are smaller 
than the bounds fj,o, fJ-i, /i2 computed in Prop. [5j Since the constants fij depend on 
the number fc, we have to rerun the program replacing fc by a real number ki larger 
than the conjectured minimum (i.e. the one obtained by running Eu3_5 on it). 
The biggest value obtained from the various exceptional points gives the Euclidean 
minimum M{K) unless they are all smaller than k. 

The situation is, however, not always as simple as in Thm. |3l In fact, look- 
ing once more at the cubic number field K with discriminant 985 and using our 
programs with k = 0.39, we can cover except for 



Ti 


[0.345,0.35] 


X 


[-0.4915,-0.49] 


X 


[-0.0185,-0.018] 


T2 


[0.0175,0.0185] 


X 


[-0.2375,0. - 2355] 


X 


[0.4725,0.475] 


T 


[0.3995, 0.4005] 


X 


[-0.201,-0.199] 


X 


[0.3995,0.4005] 


T3 


[0.4725,0.473] 


X 


[-0.146,-0.1445] 


X 


[0.2905,0.2915] 


T4 


[0.2905,0.2915] 


X 


[0.217,0.219] 


X 


[-0.437,-0.436] 


T5 


[0.436,0.437] 


X 


[0.3265,0.3285] 


X 


[0.345,0.346] 



Applying Thm. [3] to T shows that ^ — -^(2 — a + 2q;^) is the only possible 
^-exceptional point in T. Letting e = a act on the Ti we find that the 'orbit' of Ti 
is {Ti, — T2, — T3, —T4, T^}; at this point we need 

Corollary 7. Let K be a number field, e G Ek o, non-torsion unit, and suppose 
that Ti, . . . , Tt are compact subsets of F with the following properties: 
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(1) there exist /?i,...,/3f G Ok such that, for all ^ G Tj, — Pj lies in a 
k-covered region of F or in Tj+i (we put Tt+i — Ti); 

(2) for all ^ G Tj there is a j ^ Ok such that e~^^ — 7 lies in a k-covered region 
of F or in T^^d), where TT{j) is an index depending only on j and not on ^; 

(3) |e|, ^1 forl<j <r + s. 

Then the only possible k-exceptional point of Ti is C — , where j3 — e^^^Pi + 
£*~^/32 + . . . + e/3t-i + Pt- Moreover, the only possible k-exceptional points of the 
sets Tj are contained in Orbe(C)- 

Proof. Suppose that ^1 £ Ti is fc-exceptional; then ^2 = £^1 ^ Pi G T2, . . ., £,t = 
£^4_i — Pt-i G Tt, and ^4+1 — e^t — Pt ^ Ti are fc-exceptionaL Observe that 
^t+i = £^t — Pt ^ £^6-1 ~ sPt-i — Pt — ■ ■ ■ — £*^i — P- From the assumptions made 
we can deduce that, for every /c-exceptional point ^1 G Ti, £*^i — /3 hes again in Ti. 
This shows that condition 1. of Thm. [3] is satisfied with e replaced by e*. 

In order to prove that condition 2. is also satisfied we use induction to find that 
there exists a 7 G Ok such that f — e^'^(,i — 7 is an element of some set Ti. From 
what we have proved already we know that there exists a uniquely determined 
7' G Ok such that e*^ - 7' G T^. But £*f - 7' = ^1 - £*7 - 7' e implies that 
£*7 + 7' = and i = I. Thus condition 2. of Thm. [3] is also satisfied, and we can 
conclude that C = /?/(£* — 1) is the only possible fc-exceptional point in Ti. This in 
turn implies that eC — 71 is the only possible fc-exceptional point in T2, etc., and 
all our claims are proved. □ 

In our example of the cubic field of discriminant d = 985 and fc = 0.39 we now 
check that condition 1. of Cor. [7]is satisfied, and we find Pi = 0, P2 = —3a -f a^, 
Ps, = -2a, Pi ^ 2a- a'^, and P5 = 3a. This gives P = 7 -\- Ala - 28a'^ and 
^i ^ P/{a^- 1)^^(19- 27a- a^). 

After having verified condition 2., Cor. [7] shows that the only possible fc-ex- 
ceptional point of Ti is ^1 . Thus the only fc-exceptional points of X in IJ Ti are 
a = ^(19 - 27a - a2), ^ = ^(-1 + 13a - 26a^), 6 = ^(-26 + 8a - IGa^), 
a = ^(-16 - 12a -I- 24a2), ^ = ^^(24 -f 18a 19a2). Using fc = 0.5 we get the 
bounds /ii = 4.9, ^2 = 2.5, ^3 = 1-5, and M{K,^i) = N{^i + a) = f^. 

Remark 11. Suppose that, in our Example, we apply Prop. [5]to ^ — -^(2 — a-|-2a^) 
with fc = 0.39; the minimal norm of the elements satisfying conditions 1. - 3. turns 
out to be 1. Nevertheless we can only conclude that M{K,^) > 0.39. In order to 
prove that M{K,£^) = 1 we have to apply Prop. [S] again, this time with a fc > 1. 
Again, the minimal norm is 1, and now we can conclude that in fact M(Jsr, ^) = 1. 

Remark 12. Collecting the uncovered subcubes Sj into the sets Ti of Cor. [7] is 
done as follows: assume that S is an uncovered cube, and that eS — P £ F. Then 
the set Ti containing S is taken to be the set of all uncovered Si 'near' to S such 
that £5"^ — P G F. By proceeding similarly with the uncovered cubes in _F \ Ti we 
eventually arrive at subsets Ti containing all uncovered subcubes. 

The programs are available from the authors. We remark that they can also be 
used to study weighted norms; details will be presented in [3]. 
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4. Some Heuristic Observations 

Consider some ^ = ^ + Ok G K/Ok', the bigger # Orb(^), the more hkely it is 
that one of the points in the orbit can be approximated sufficiently well by some 
7] e Ok- In fact, if (a,/3) = (1) and if #Orb(0 is maximal (i.e. #Orb(C) = 
{Ok ■ POk)^), then clearly M{K,^) = \/N(i. Euclidean minima thus tend to be 
attained at points ^ with small orbits. If K has unit rank 1, then there are many 
points with small orbit: just take any for a S Ok and e G Ek a fundamental 
unit. If the unit rank is > 2, however, such points are hard to find, because there 
is no guarantee that a/{ei — 1) has a small orbit with respect to the action of a 
second unit £2. 

There is one exception, however: suppose that there is a principal prime ideal 
p — (tt) which is completely ramified in iC/Q. Since p has degree 1, for every 
e G Ek there is an integer a G Z such that e = a mod p. Taking the norm 
gives ±1 = Nk/qs = Nk/qcl ~ a" mod p, where n — {K : and this in turn 
implies that e" = a" = ±1 mod p. Therefore the unit group generates at most 2n 
different residue classes mod p, hence # Orb(^) < 2n for any ^ of type ^ = ^ + Ok- 
Therefore such ^ have comparatively small orbit and a good chance of producing a 
large minimum. In fact, almost all known Euclidean minima of normal cubic fields 
are attained at such points. 

Another question we would like to discuss is the following: can we expect that 
our list of norm-Euclidean complex cubic number fields is complete? Let us see 
what is happening in the real quadratic case. There we know (cf. f8|) that (in the 
following, d = disc K denotes the discriminant of K) 

^ < M{K) < ^ 



16 + 6v/6 ~ ' ' - 4 ■ 

This allows us to define the Davenport constant D — sup AI{K)/^/d for real qua- 
dratic fields. The example K = Q(\/l3), M{K) = 1/3 shows that D > l/SvTS. 
If we assume that this is a good approximation for D, then there should be no 
norm-Euclidean number fields with discriminants > D^^ = 9 • 13 = 117; in fact, 
the maximal discriminant of a norm-Euclidean number field is d = 76. 

If we try to do the same with complex cubic case then the first problem is that 
the exponent 1/2 of the discriminant in the lower bound in 



I''! < MiK) < M^'" 



420 - ' ' - 

is not known to be best possible. If it is, then we can define a Davenport constant 
D — sup M{K)/\/d for complex cubic fields as well. The example d — —244, 
where M{K) = 1/2, shows that D > 1/2^244, and if this bound is good, then 
there should be no norm-Euclidean number fields with \d\ > 976. The example 
d = —503 suggests that D is somewhat smaller, but in any case we don't expect to 
find norm-Euclidean fields with \d\ > 1500. Basically the same conclusions (with 
better bounds) hold if the correct exponent of \d\ in the lower bound is 2/3. 

In the case of totally real cubic fields there is no known (nontrivial) lower bound 
for M{K) at all (of course M{K) > |). If one could show M{K) > cVd for some 
c, then the above heuristics show that one has to compute M{K) at least for fields 
with discriminants up to 25. 000 (in fact Godwin and Smith [6] have shown that 
the normal cubic field with discriminant d = 157^ = 24.649 is norm-Euclidean); 
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our current data are therefore insufficient for deciding whether such a lower bound 
might exist or not. 

5. A CONJECTURE 

We would like to conclude our paper with a conjectur^ concerning M{K) for 
certain pure cubic fields K: 

Conjecture 1. Let to = + 1 be a squarefree integer, and assume that I is even; 
put a = VI, K ^ Q(a), and ^ = i(l + a + a^). Then 

M(K) . MiK, = l¥'''l ' K ^ ^ 

^ ' ^ ^' [ ^(18^4 _ _^ 30^2 _^ -32), if^ = 0mod4. 

It is easy to see that M{K,£^) has at most the value given above; in fact, if 
i = 2 mod 4, then N{l£^ + \ + \£a- ^a^) = ^(18^"* - + 12^2 + 12^), and if 
£ = mod 4, then N {\P + \i+\ + {\l-\)a-\a^) = ^{lU^ -U^ +2,0^ +2Al-'i2) . 
Numerical computations show that the conjecture is true for £ = 4 and H. ~ 10. 
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These Tables contain the known Euclidean minima for cubic number fields of 
small discriminant. The fields are ordered by |disc-ftr|; fields with equal discrimi- 
nant are ordered as in the number field tables at Bordeaux 
(file : //megrez . ceremab .u-bordeaux . f r/pub/numberf ields/). 
The letter N indicates that the field has class number 1 but is not norm-Euclidean, 
and H that it has class number > 1. Moreover, E means that the field is norm- 
Euclidean; if no Minimum is given, we succeeded in covering the fundamental do- 
main with k — 0.99. CPU-times ranged from a few minutes for fields with small 
discriminants to several hours; by far the hardest nut to crack was AiscK — 10661, 
which took several days. 

The only fields with disc if < 11,000 whose Euclidean nature is currently not 
known are those with discriminants 10929 and 10941. We also remark that among 
the four fields which were shown to be Euclidean in [6], those with discriminants 
11881, 16129 and 24649 are beyond the limits of our tables. 

(*) The Euclidean minimum M{K) for the field with disc if = —680 is 

81956632 
- 81182612- 

(§) The Euclidean minimum M{K) for the field with disc if — —728 is 

7483645229 



M{K) 



8158377554' 



